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INTRODUCTION 


It  is  well  known  that  the  closed-form  elastoplastic  solutions  of  the 

generalized  plane-strain  problems  can  only  be  obtained  by  Hencky’s  total 

1,2 

stress-strain  relations  or  Tresca’s  yield  criterion.  Based  on  the  detailed 

studies  reported  in  reference  3,  more  realistic  material  model  for  thick- 

walled  tube  under  very  high  pressure  operation  is  an  elastic-plastic  material 

which  obeys  the  Mises'  yield  criterion  and  the  Prandtl-Reuss  incremental 

stress-strain  relations.  And,  in  such  situations,  one  has  to  rely  on 

4~  6 

numerical  methods.  Both  the  finite-difference  method  and  the  finite- 
7,8 

element  method  have  been  used  to  solve  the  elastoplastic  problem  considered 

here.  The  finite-element  method  is  more  powerful  and  can  be  used  to  solve 

8“10 

more  general  elastoplastic  problems.  Since  the  displacement  function  is 

assumed  and  the  programming  is  complicated,  the  accuracy  of  the  finite- 
element  approach  has  to  be  verified.  This  is  usually  done  by  comparison  with 
more  rigorous  solutions  to  simpler  problems.  For  the  problem  considered  here, 
rigorous  solutions  based  on  the  finite-difference  method  were  obtained  for 
ideally-plastic  materials  in  references  4  and  5  and  for  strain-hardening 
materials  in  reference  6. 

In  the  present  paper,  a  new  finite-difference  approach  is  developed  for 
solving  the  generalized  plane-strain  problems  of  thick-walled  cylinders 
subjected  to  internal  pressure,  external  pressure,  or  end  force  beyond  the 
elastic  limit.  An  incremental  approach  is  used  and  no  iteration  is  needed  for 
each  increment.  The  numerical  scheme  is  stable  for  ideally-plastic  as  well  as 

References  are  listed  at  the  end  of  this  report. 
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strain  hardening  materials.  The  approach  is  simpler  than  others  yet  quite 
general  and  accurate.  The  desired  accuracy  can  be  achieved  by  reducing  the 
grid  size  and  load  increments. 

BASIC  EQUATIONS 

Assuming  small  strain  and  no  body  forces  in  the  axisymmetric  state  of 
generalized  plane  strain,  the  radial  and  tangential  stresses,  ar  and  Oq,  must 
satisfy  the  equilibrium  equation, 

r( 3ar/9  r)  =  a0  -  ar  ;  (!) 

and  the  corresponding  strains,  er  and  e0,  are  given  in  terms  of  the  radial 
displacement,  u,  by 

e  r  =  3u/3r  ,  e0  =  u/r  .  (2) 

It  follows  that  the  strains  must  satisfy  the  equation  of  compatibility 

r(9e0/9r)  =  er  -  e0  .  (3) 

Whereas  the  differential  equations  (1),  (2),  and  (3)  hold  throughout  the  tube 

regardless  of  the  material  properties,  the  constitution  equations  assume 

various  forms  according  to  the  adopted  form  of  yield  function,  hardening  rule, 

total  or  incremental  theory  of  plasticity.  In  the  present  paper,  the  material 

is  assumed  to  be  elastic-plastic,  obeying  the  Mises'  yield  criterion,  the 

Prandtl-Reuss  flow  theory  and  the  isotropic  hardening  law.  The  complete 

5 

stress-strain  relations  are: 

dEi'  =  dai'/2G  +  (3/2)ai'da/(oH')  (4) 

da  >  0  for  i  =  r,Q,z 

d£m  =  E-1(l-2v)dam  (5) 


2 


where  E,  v  Young* s  modulus,  Poisson* s  respectively, 

2 G  =  E/(l+v) 

cm  =  (£r+£0+£z)/3  >  ei'  =  ei  ~ 


"m  > 


am  =  (or+a0+az )/3  »  °i'  "  ai  -  am  » 

a  =  ( 1//2) [ (ar-O0)2  +  (ae-az)2  +  (az-ar)2 ] !/2  >  a0  ,  (6) 

and  oQ  is  the  yield  stress  in  simple  tension  or  compression.  For  a 
strain-hardening  material,  H*  is  the  slope  of  the  effective  stress/plastic 
strain  curve 


o  =  H(JdeP)  .  (7) 

For  an  ideally-plastic  material  (H*  =  0),  the  quantity  ( 3/2)da/ (aH* )  is 

replaced  by  dX  ,  a  positive  factor  of  proportionality.  When  a  <  oQ  or  da  <  0, 

the  state  of  stress  is  elastic  and  the  second  term  in  equation  (4)  disappears. 

9 

Following  Yamada  et  al,  equations  (4)  and  (5)  can  be  rewritten  in  an 
incremental  form 


dOf  =  d-j_jd£j  for  i,j  =  r,0,z 

and 

dt j/2G  =  v/(l-2v)  +  -  Oj/OjVS  ,  (8) 

where 

2  1 

S  =  -  (1  +  -  HT/G)a2  ,  H'/E  =  a/(l-a)  ,  (9) 

aE  is  the  slope  of  the  effective  stress-strain  curve,  and  is  the  Kronecker 

delta. 

This  form  was  used  in  the  finite-element  formulation  for  solving 

7 

elastic-plastic  thick-walled  tube  problems.  In  the  following  section,  the 
incremental  stress-strain  matrix  will  be  used  in  the  finite-difference 
formulation. 
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FINITE-DIFFERENCE  FORMULATION 


Consider  an  open-end  or  closed-end  thick-walled  cylinder  of  inner  radius 
a  and  external  radius  b.  The  tube  is  subjected  to  inner  pressure  p,  external 
pressure  q,  and  end  force  f.  The  elastic  solution  for  this  problem  is  well- 
known  and  the  pressure  p*,  q*,  or  f*  required  to  cause  initial  yielding  can  be 
determined  by  using  the  Mises*  yield  criterion.  For  loading  beyond  the 
elastic  limit,  an  incremental  approach  of  the  finite-difference  formulation  is 
used.  The  analysis  starts  with  the  applied  pressure  p,  q,  or  f  and  the  load¬ 
ing  path  is  divided  into  m  increments  with 

Ap  =  (p-p*)/m  ,  Aq  -  (q-q*)/m  ,  Af  =  (f-f*)/m  .  *  (10) 

The  cross  section  of  the  tube  is  divided  into  n  rings  with 

rl“a>r2> • • • >rk”P > • • • »rn+l*b  >  (11) 

where  p  is  the  radius  of  the  elastic-plastic  interface.  At  the  beginning  of 
each  increment  of  loading,  the  distribution  of  displacements,  strains  and 
^stresses  is  assumed  to  be  known  and  we  want  to  determine  Au,  Aer,  Agq,  Aez, 
Aor,  Aoq,  A oz  at  all  grid  points.  Since  the  incremental  stresses  are  related 
to  the  incremental  strains  by  the  incremental  form  (Eq.  (8))  and  Au  »  rAeQ, 
there  exists  only  three  unknowns  at  each  station  that  have  to  be  determined 
for  each  increment  of  loading.  Accounting  for  the  fact  that  the  axial  strain 
ez  is  independent  of  r,  the  unknown  variables  in  the  present  formulation  are 
(Agq)^,  (Aer)^,  for  i  =  1 ,2, . . .n,n+l ,  and  Aez. 

The  equation  of  equilibrium  (1)  and  the  equation  of  compatibility  (3)  are 
valid  for  both  the  elastic  and  the  plastic  regions  of  a  thick-walled  tube. 

The  finite-difference  forms  of  these  two  equations  at  i  =  l,...,n  are  given  in 
reference  6  by 


4 


C r i_|_ i~ 2ri) (Aar) i  -  (ri+1-ri)(Aa0)i  +  ri(Aar)i+1 

=  (ri+i-riHag-Oj.^  -  ri[(ar)1+1  -  (ar)i]  (12) 

for  the  equation  of  equilibrium,  and 

( r i+l-2r i)  ( A£q )i  -  ( ri+l_r i)  ( ^er^  i  ri(Ae0)i+l 

=  ( r i+i—r i )  ( £ g q ) i  —  r-^ [  ( £ q ) ^ q ) i ]  (13) 

for  the  equation  of  compatibility.  With  the  aid  of  the  incremental 
s tress-strain  relations  (Eq.  (8)),  equation  (12)  can  be  rewritten  as 

[ Cri+l-2ri)(d12)i  +  (-ri+l+ri)(d22)il(Aee)i 
+  [(ri+i~2ri)(d11)i  +  (_ri+1+ri) (d2i)  j.1  ( AEr) i 

+  ri(d12)i+i(A£e)i+l  +  ri(dll)i+l(AEr)i+l 

+  [(ri+1-2ri)(d13)i  +  (-ri+1+ri)(d23)i  +  ^i(d13) i+1 ] Aez 

=  ( ^i+i— ^i) ( ®0— 1 ®  r) i  —  —  *  (1^) 

The  boundary  conditions  for  the  problem  are 

Aar(a,t)  «  -Ap  ,  Aar(b,t)  =  -Aq  , 
n 

n  l  [ri(Aaz)i  +  ri+1(Aaz)i+1]  (r^-r^  =  inra2Ap  +  Af  ,  (15) 

i=l 

where  y  is  0  for  open-end  tubes  and  1,  for  closed-end  tubes.  Using  the 
incremental  relations  (Eq.  (8)),  we  rewrite  equation  (15)  as 

(d12h(Ae0h  +  (dll^l^Aer)l  +  ^d13^1Aez  =  ~AP  »  (16) 

(d12)n+l(Aee)n+l  +  ^dl  Pn+l  (  Aer)n+1  +  (dl3)n+lA£z  =  _A<1  >  U7) 

and 

l  (ri+l“ri){  rit(d23^i^Ae0^i  +  (d13^ i< A£r> i.3  +  ri+l  ^  (d23)  i+1  ^^i+l 
i=l 

n 

+  (dl3) i+l (Aer^ i+1 1 }  +  I  (ri+l~ri)  tri(d33) i  +  ri+l^d33)l+l] Aez 

=  pa2Ap  +  Af/ir  .  (18) 
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Now  we  can  form  a  system  of  2n+3  equations  for  solving  2n+3  unknowns,  (Aeq)^, 
(Aer)i,  at  i  =  1 ,2, . . . ,n,n+l  and  Aez.  Equations  (16),  (17),  and  (18)  are 
taken  as  the  first  and  last  two  equations,  respectively,  and  the  other  2n 
equations  are  set  up  at  i  =  1,2,... ,n  using  (13)  and  (14).  The  final  system 
is  an  unsymmetric  matrix  of  arrow  type  with  the  nonzero  terms  appearing  in  the 
last  row  and  column  and  others  clustered  about  the  main  diagonal,  two  below 
and  one  above. 

In  the  computer  program  which  was  developed,  the  dimensionless  quantities 
r/a,  Eer/o0,  Eeq /o0,  Eez/o0,  ar/o0,  og/a0,  az/o0,  p/aQ,  q/a0,  f/(7ra2a0)  were 
used  in  the  formulation  and  the  Gaussian  elimination  method  was  used  to  solve 
these  equations.  All  calculations  were  carried  out  on  IBM  360/Model  44  with 
double  precision' to  reduce  round-off  errors. 

.NUMERICAL  RESULTS  AND  DISCUSSIONS 

The  generalized  plane-strain  problems  of  thick-walled  cylinders  subjected 
to  internal  pressure  p  beyond  the  elastic  limit  were  solved.  The  elastic- 
perfectly-plastic  as  well  as  strain-hardening  materials  were  considered  for 
open-end  or  closed-end  conditions.  The  numerical  results  were  based  on  the 
following  parameters:  b/a  =  2,  v  =  0.3,  a  =  0.05,  p  =  0  or  1.  Various  values 
of  m  and  n  were  used  to  test  the  convergence  of  the  numerical  results.  The 
incremental  loadings  were  applied  until  the  fully  plastic  state  was  reached. 
The  value  for  p  corresponding  to  this  final  state  was  denoted  by  p**.  It  was 
found  that  the  results  of  these  values  for  all  four  cases  converge  by 
increasing  m  and/or  n.  For  n  =  100,  Ap/a0  =  0.0004,  we  have  p**/oQ  =  0.7802 
(a  =  0,  p  =  0  as  case  1);  0.8038  (a  =  0,  p  =  1  as  case  2);  0.8222  (a  =  0.05, 
p  =  0  as  case  3);  0.8618  (a  =  0.05,  p  =  1  as  case  4).  Additional  results  are 
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shown  in  Figures  1  to  5.  Figure  1  shows  the  bore  radial  displacements  as 
functions  of  internal  pressure  for  cases  1,  2,  and  4.  Figure  2  shows  the 
relations  between  internal  pressure  and  elastic-plastic  boundary  for  cases  1, 
2,  and  4.  The  effects  of  end  conditions  and  strain  hardening  can  also  be  seen 
in  these  two  figures.  The  distributions  of  radial  and  tangential  stresses  for 
p /a  —  1.0,  1.2,  1*4,  1.6,  and  1.8  are  shown  in  Figure  3  for  case  1  and  in 
Figure  4  for  case  2.  Finally  the  distributions  of  axial  stress  for  p /a  =  1.0, 
1.4,  and  1.8  are  shown  in  Figure  5  for  cases  1  and  2.  The  effect  of  end 
conditions  and  elastic-plastic  boundary  on  the  axial  stress  is  quite 
significant. 

The  present  approach  determines  Aez  directly  for  each  step  of  incremental 
loadings  whereas  in  reference  6,  many  iterations  were  needed  because  a  value 
of  Ac z  was  assumed.  In  addition,  the  computer  storage  needed  in  this  approach 
was  only  35%  of  that  in  reference  6,  and  much  larger  n  can  be  used  to  yield 
better  results.  The  present  approach  is  simpler  yet  more  general  than  the 

4,5 

other  finite-difference  approaches  because  both  ideally-plastic  *  and  strain- 
6 

hardening  materials  can  be  considered  in  a  unified  manner.  Furthermore,  very 

accurate  numerical  results  can  be  obtained  and  used  to  verify  the  accuracy  of 

7,8,10 

the  finite-element  programs. 
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Figure  1.  Radial  displacement  at  the  bore  as  a  function  of  internal  pressure  for  cases 
1,  2,  and  4. 
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Figure  2.  Elastic-plastic  boundary  as  a  function  of  internal 
pressure  for  cases' 1,  2,  and  4. 
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Figure  3.  Distributions  of  radial  and  tangential  stresses  for 
case  1  (a  =  0,  y  =  0,  b/a  =  2). 
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Figure  4.  Distributions  of  radial  and  tangential  stresses  for 
case  2  (a  =  0,  y  =  1,  b/a  =  2) . 
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